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1 Introduction 



The QCD description structure of the final state in high energy processes is one of the 
most important subjects in strong interaction physics. Of special interest is the analysis 
of the final state formation at high energies, where one should resum the contributions 
logarithmic in energy (Bjorken x). This problem has drawn a considerable interest and 
was studied in a number of recent publications 0- 

In the approach traditionally applied for describing the final state in the high energy 
hadron collisions, that of coUinear factorization, the final state is produced with the zero 
total transverse momentum. In particular for two-particle final state this means that these 
particles leave the collision point in the opposite directions and having equal absolute 
values of their transverse momenta. We will refer to such a configuration as symmetric 
and will consider any deviation from it as asymmetry. The value of the asymmetry can 
be used as the measure of the departure from the collinear factorization showing limit 
of its applicability. To quantify the effects related to the nonvanishing total transverse 
momentum of the final state one has to use an approach in which the transverse momenta 
of incoming partons are not neglected. Let us note, that in a number of recent experiments 
on diphoton 0, 7r° and direct photon |^ production a substantial discrepancy between 
the data and predictions of coUinearly factorized NLO QCD was observed. Taking into 
account the intrinsic transverse momenta of order of 1 — 2 GeV substantially improves 
this situation. The results of 0| also show sufficiently broad angular distribution in 7r° 
production. 

The paper is organized as follows. In the next section the high energy factorization will 
be reviewed and compared with the collinear factorization one. In section |^ numerical 
results for the shape of the final state containing two particles in the central rapidity 
region are presented. The derivation of the expressions for the cross sections used in our 
calculations can be found in the Appendices. In section ^ we summarize the results and 
present our conclusions. 

2 Collinear and high energy factorization 

Collinear factorization |^ is a method of describing the strong interaction processes by 
factorizing the contribution to physical cross sections into the product of partons distri- 
butions fa{x,k'^) parametrizing both the nonperturbative information about the hadron 
and perturbative evolution starting from some specific initial condition and perturbative 
cross sections corresponding to the scattering of the parton fiuxes. In this approach the 
prehistory of colliding partons is entirely determined by structure functions while partons 
themselves are supposed to be on-shell particles. The cross section for two-particle (jet) 
production to the lowest perturbative (Born) order is 




da 



da, 
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where the sum over parton types a, 6 is assumed and 0:1,2 = /c±(e^^^ + e^^^)/\fS, where 
\/S is an invariant colhsion energy. 
In this approach it is assumed that 

Aqcd < fc^ ~ xy/S ~ v^. (2) 

indicating that one can apply perturbative QCD and that there is only one big logarithm 
to take into account, that is hi{k'^ / Resummation of the powers of this logarithm 
(i.e. ~ a^WiJi^ / ls?Q(jj-))) leads to structure functions depending on k"^ and is performed 
by DGLAP evolution equation 0. 

The situation changes at large S when one reaches the kinematic region 

Aqcd < ~ x\fS < (3) 

Under these conditions another big logarithm ln(l/x) exists. Resummation of powers of 
this logarithm can become more important than the one of ln(A;^/AQ(^^). The resum- 
mation of the leading energy logarithms for the structure function is described by BFKL 
equation |^. 

In this kinematical region the transverse momenta of the incoming parton fluxes can no 
longer be neglected. To take them into account a new approach called " factorization" 
was proposed in |^ . Extensive description of the method and various applications can be 
found in 0. Let us note that this method was de facto used earlier in |T0[| . 

The method of k^ factorization is based on considering "partons" with nonzero trans- 
verse momentum and being, in contrast to the collinear factorization case, off-mass shell 
particles. Calculation of physical cross-sections requires a generalization of the conven- 
tional picture of parton flux as described by usual parton structure functions. To this aim 
an unintegrated gluon distribution (j){x, qj_, k) is introduced: 

-^(f){x,q±,k), (4) 
1± 

where is proportional to the probability of finding the incident parton with longi- 
tudinal momentum component xpa {pa is a momentum of initial particle) and transverse 
momentum q±. It is important to stress that unintegrated structure function depends not 
only upon the momentum of particle to which it corresponds, q, but also on the global 
off-mass-shellness of the process, fc^. This is in fact a general property of distributions in 
quantum theory - the probability of a particular event depends not only of its parameters, 
but also on the global characteristics of the event ensemble. 

The particular interrelation between the unintegrated and integrated structure func- 
tions depends on the evolution equation that the integrated distribution solves. For the 
DGLAP evolution ^ the unintegrated structure function can be expressed through the 
integrated one as follows: p3l 



(x, q±, k) 



q;„ t dz 



X 



J -^P,,{z)xg(^^,qlj\T,{q^,k) (5) 
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where Tg is the Sudakov form factor for gluon 



Tg{q±, k) = exp 



27r Pi 



i-=g,<},<i 



(6) 



(for Pgg in the last expression one should substitute 



99 \ 



ZPgg{Z 



In a double logarithmic approximation Eqs. 



in account of 
coincide with the 



gluons identity). 
DDT formula 

For BFKL evolution 0] the correspondence between integrated and unintegrated struc- 
ture functions takes a simpler form 



{x,q±) = (p{x,q±,q±) 



dxg{x,ql) 
dlnqi 



(7) 



This results from the absence of strong q± ordering in BFKL ladder. 

From the practical point of view, Eqs. (0-^ differ substantially from Eq. (|^) only for 
q± k implying the presence of a big logarithmic contribution oc In/c^/g^. Therefore we 
can in fact use Eqs. (|1-|D for both types of QCD evolution. 

The calculation of a cross section with virtual colliding particles presents a highly non- 
trivial problem. The first difficulty lies in the correct summation over polarization states 
of the incident off-shell "partons". Another problem is the presence of bremsstrahlung 
from initial and final states of colliding particles. The derivation of the relevant cross 
sections is outlined in Appendix A. 

More specifically, the transverse (high energy) factorization is based on considering the 
2 ^ n + 2 process cross section with a large rapidity gaps between the two final particles 
providing a full rapidity interval for the process under consideration which are almost 
collinear to the incident ones and n particles emitted into the central rapidity region (so 
called quasi-multy-Regge kinematics or QMRK 0)- Note that only the amplitude 
with on-shell and physically polarized in- and outgoing particles has precise gauge invari- 
ant meaning and only such expressions can be used in calculating the cross section. Large 
rapidity gaps provide a separation between quantities related to the incident particles and 
those describing the hard cross section of parton production in the central region. 

Let us for example consider the process gg —>■ ggg in the limit of high energy 




99 999 
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with gi,2 = Pa,b - p'a,b- 

This equation demonstrates the above mentioned factorization of the cross section. 
The first, second and third factors under the integral correspond to Pa p'ai Qi sphtting, 
Qi,(l2 ^ k "scattering" and pb — > p'2, q2 sphtting respectively. 

The factors related to the splitting of the incident particles should further be converted 
to structure functions. One can do it in two steps. The first step, transformation to the 



form factors, is quite straightforward [|13[. In the second step we have to account for 
corrections due to radiation along the directions of incident particles and replace form 
factors by unintegrated structure functions (p{x, q±, k) with x determined by kinematics. 
The cross sections for producing n = 0, 1, 2 particles in the central region read 

(27^)^ [ .2 v{.x,qi_)ip{x,qA_) 



- / d q±—o -2 X = qJS; (9) 

I J g_|_ q^ 

/( ' 
d^qi±d'^q2j_- 

(10) 



d'^k±dy J ^"^"^ ^^"^ qlj_ dkj_ g|_|_ 
dai ANctts 5(2) + g2± - k±) 



dkl kl 
xifi = k±_e^ /\/~S, X2fi = k±e~^ /\/S] 




99 9999 



dao 



d'^ki±d'^k2±dyidy2 



^ J d'^qi±d'^q2± 



(p{xi,qi±,k±) da2 ^p{x2,q2±,k±) 



ddo 



d'^ki±d'^k2± 
k± = ki± + k2± 
2Nlal 5(2) (g^^ + - ki^ - k2± 



A, 



d^ki^d^k2± - 1 

xi = (ki^e^' + k2±ey')/y/S, X2 = {k^^e-^^ + k2^e-y^) / ^/S . 



'Ill 
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The expression for n = corresponds to the total cross section in the two gluon exchange 
(Low-Nussinov) approximation. Let us note that the kinematics of Eq. (|^) is not quite 
the one in Eq. @. Having no kj_ in this process we have to restrict directly qj_. In Eq. (P) 
q± = \/xS while it is q± ~ x^/^ that is given by (^. 

The case of one gluon production in the central region {n = 1) was studied in a number 
of publications, e.g. |T^, [p!4| . 



A part of the analytical expression for A (for gg — > gg subprocess) has recently 
been published in [0, |T5|. A similar quantity for gg qq subprocess can be found in 



T^, |T6|. The derivation of the explicit expressions is outlined in the Appendices A and 
B. Let us note that the formula for gg qq cross section from |TB[ coincides with the 
analogous formula in (in the limit of massless quarks) obtained by direct application 
of k± factorization. 

For practical applications of the equations Eq. (P-|TT|) we should understand the overall 
normalization of the off-shell cross sections to the usual on-shell one arising in collinear 
factorized formalism. The normalization can be deduced from considering the small 
limit in Eq. (0). Physically, we have to return to Eq. (§) and choose the structure 
functions being extremely narrow functions of q\_ and integrate out (Pqii_ and (Pq2±- There 
is no contradiction with Eq. (H) because q± <^ k±. Certainly, we must put gi_L,g2± = 
in (3"2 Eqs. (|I^, [B.2| , |B.3D . As it is follows from Appendix A, A is proportional to 
so this limit for (3"2 is correct. After this substitution and integration over k2± 
Eq. (|ip appears with the structure functions from Eq. and with the correct da/dt. 
Consequently, normalization of Eqs. (|9|-|TT|) is also correct. Let us emphasize that the 
integration over d'^qi± and d'^q2± includes averaging over angular orientations of qi± and 
q2± in the transverse plane when arriving to the final expression for o"2. This averaging 
is similar to the averaging over initial gluon polarizations in getting the usual expression 
for the cross section of the 2 — > 2 elastic scattering. 



3 Angular and momentum asymmetry of parton pro- 
duction at central rapidity 

The most interesting predictions of the high energy factorization are of course those going 
beyond the collinear factorization framework. The simplest situation corresponds to the 
production of a single particle in central rapidity region described in the first order in as 
Eq. (0). This process was first studied in |jl3| and later in many publications, see e.g. 



0- 

In this paper we concentrate our attention on another important feature of the high 
energy factorization. From Eq. (|TT|) it is obvious that, in contrast with the collinear fac- 
torization case, the two partons produced in the central rapidity region are not necessarily 
going into opposite directions (back-to-back). Another important issue is that the abso- 
lute values of the momenta of outgoing partons are not necessarily equal. We will refer 
to these properties as angular and momentum asymmetries respectively. 

The same characteristics of the two-particle production were studied in [|T^ in the case. 
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where the produced particles are separated by the large rapidity gap. The angular and 
momentum asymmetries considered in are due to the presence of BFKL ladder filling 



the gap between the produced particles. Let us stress, that in this paper we study the 
two-particle decorrelation in a relatively small central rapidity interval. The rapidity gap, 
discussed in , is thus absent and decorrelation effects are exclusively due to off-shellness 



of the two generalized parton fluxes merging in a particle-producing vertex. 

Our approach is similar to that in [^, where a multy-regge limit for A was considered, 
A = (li±(l2±/^i±^2± (cf- Appendix B). In this form A has no intrinsic angular and 
momentum correlations (see below). Moreover, the contributions from A and the structure 
functions to CT2 (see Eq. (0)) factorize and (T2 depends on structure functions only via 



gluon-gluon luminosity function p^. The latter contains practically all information about 
particle correlations. 

A large number of independent variables in Eq. ([TT|) makes it difficult to visualize the 
pattern of particles emission. Therefore we will use the less differential ones, in particular 
the angular asymmetry 



dA(f) dyidy2 



^ f Jr, , dki^dki^, (12) 

J d^k, I d^ko idv^dvo 



fco 



d'^ki±d'^k2±dyidy: 



where A(j) G [0, vr] is the angle between particles. Integration over ki± and k2± requires 
introducing the infrared cutoff limit /cq. The shape of asymmetry is sensitive to the chosen 
value of ko- In the normalization of Eq. (|T2|) and in forthcoming Eq. (|T5D we took into 
account particles identity (even for quarks we do not distinguish between q and q). 

In order to perform actual calculations one needs to choose the unintegrated structure 
functions entering Eqs. (|9|-p!ID. In principle they should be chosen as solutions to the 
NLO BFKL equation [|18] or its nonlinear generalizations [T^ (also taking to account 
form factors arising from non-compensating real and virtual corrections). Leaving this 
analysis for the future, let us note, that the minimal requirement surely has to be that 
whatever expression for the unintegrated distribution is used, it should not contradict the 
information from deep inelastic scattering, i.e. that on integrated structure function. This 
fixes the distribution in the limit q± = k±, where 4>{x,q±,q±) = 4>{x,q±) (cf. Eq. (|^)). 
In our numerical calculations we have used the unintegrated gluon structure functions 
corresponding to integrated CTEQ5M |^T| and AKMS p2| fitted as in |jl^. For q±_ < k± 



the unintegrated structure function (j){x, q±, k^) was determined according to Eqs. (||), (||). 
In the final distributions we sum over gluon and quark {uf = 5) pair contributions (in 
fact, the contribution due to quark pair production is less than 3% of that due to gluons). 

In Fig. |l| we show differential cross section (|12D with yi^2 = yoT ^Z//2 for yo = and 
Ay G [0, 2] and with ko = 2GeV. From now on figures marked by a, b are calculated 
with CTEQ structure functions, marked by c, d are calculated with AKMS structure 
functions; figures with labels a, c correspond to = l.STeV, while those with labels b, 
d correspond to -y/S* = 14TeV. 

The most striking feature of this cross section that differs it from the collinear factor- 
ized one is the appearance of the collinear singularity at Ay, A0 — > 0. This singularity is 
just a well-known s-chanel one and is in turn a reflection of the presence of one particle 



7 



production process Eq. {^^. Moreover, the behavior of the two-particle production ( [Tl| ) 
near the above- described coUinear singularity can be presented in the following decom- 
posed form. Introducing = Ay^ + A0^, = ki± + k2±, z = ki/k (since ki and ^2 are 
collinear this definition is unambiguous) and taking the limit r — > 0, we get 

l"'-(P„(.) + 2,vP„(.)), 



(Pk^dyQcirdz (Pk±dyQ r vr 



with one particle production cross section as defined in (|iy). Pgg and Pgg are the standard 
Altarelli-Parisi kernels [^] (the factor of 2 before Pgg is due to identical treatment of quarks 
and antiquarks in our approach). From Fig. |l] we can conclude, that the dependence of the 
asymmetry on the structure function is quite pronounced, see corresponding upper and 
lower panels. For AKMS structure functions the angular distribution is noticeably wider 
then for CTEQ ones. This is the reflection of the fact that AKMS structure function is 
broader then CTEQ one if compared at some fixed value of x. The dependence of angular 
asymmetry on c.m.s. collision energy is, as seen from Fig. |1|, relatively weak. 

In Fig. 0we take a closer look at the angular asymmetry by plotting the two-dimensional 
cross sections of the plots in Fig. |I] for fixed values of Ay, i.e. study the dependence of 
(|1^) on for different values of Ay where 

/ » ,N da-? 
p{A(t>) 



dA(j) dyidy2 
is normalized according to 

^^P(A^) = 1 (13) 

71 

for each value of Ay. A substantial deviation of p from the back-to-back shape, for which 
p ~ 6{A(j) — tt), is obvious. From Fig. |^ we see that the asymmetry grows with increasing 
collision energy, although the growth is not pronounced. Using different structure func- 
tions affects the shape of the asymmetry distribution much stronger than changing the 
collision energy. 

When the cutoff ko increases, the kinematic interval Eq. (^) narrows and the two 
outgoing particles tend to appear back-to-back. This dependence is clear from Fig. ^ 
where p(A0) is shown for Ay = 1 and several values of fco- 

We see that when k± ^ Qchar, where qchar is a characteristic transverse momentum 
carried by the parton fluxes, jets are predominantley produced at small 50 = vr — A0. 
This regime was studied in DLA in [^. According to ||25[ the (normalized) azimuthal 



asymmetry in 50 <^ 1 domain is given by 

with Tg{6(l)) = Tg{6(f) ko, ko) (cf. Eq. (H)). From Eqs. (|]) and ([I^) it is obvious that 
the rapid growth of p with diminishing 50 stops at 50* ~ expl—n/Gas)- Numerically, for 
ttg = 0.18 (this value of strong couphng was used in CTEQ5M calculation for k = 20GeV), 
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(50* ~ 5 ■ 10~^, in agreement with Fig. |^. Note that contrary to DLA in our approach the 
value of p for 5(f) < 5(f)* is not decreasing. The reason is that in our case there is no scahng 
leading to Eq. (|T4|). In Fig. |^ the effect of saturation is manifest not for p{5(f)) itself, but 
for its derivative, i.e. for the corresponding characteristic angle. 

Turning now to the analysis of the momentum asymmetry let us consider the cross 
section integrated over the angular variables, 

2hk, I ^ d(f),d(f),. (15) 



dki±dk2±dyidy2 J d'^ki±d'^k2±dyidy2 

For yi^2 = ±0.5 and A;i.2± G [2, 20] GeV this cross section is plotted in Fig. ^. From this 
figure we see that, in a characteristic event, ki± and k2± are not equal. The dominant 
trend can be described as a decrease of the cross section with growing k±. However, 
for the values of fci,2± of order of, or higher than 10 GeV, the correlation pattern is no 
longer powerlike. Finally, let us note that the momentum asymmetry shows a significant 
dependence on the choice of the structure function and on c.m.s. collision energy. 



4 Conclusions 

Angular and momentum asymmetry is a characteristic feature of particle production in 
semihard kinematic region. When the momenta of produced particles are of the same order 
as the characteristic transverse momenta carried by the generalized off-shell partonic fluxes 
described by the unintegrated structure function, the angular distribution of particles 
shows significant deviations from the conventional back-to-back picture. The asymmetry 
dies away when the transverse momenta of produced particles are, correspondingly, much 
larger than the intrinsic ones. The momentum asymmetry pattern is somewhat different 
due to rapid growth of the cross section with decreasing momenta. Nevertheless some 
reflection of momentum balance in the case of relatively high momenta is shows itself 
through local cross section deviation from power-like regime around the point where the 
transverse momenta of produced particles are equal. 
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Appendix A Particle production in QMRK 



As mentioned before the QMRK regime is that in which the incident particles scatter 
at parametrically small angle producing particle(s) in the central rapidity region. The 
leading contribution to the scattering amplitude in this kinematics has the form (in the 
Feynman gauge) 

A2^n+2 = /r-4p:^^ivf;;;>rlr-, (a.i) 

HI H2 

where F is a (helicity conserving) vertex and i stands for the adjoint representation index. 
The incident particles have initial momenta pa and ph and the final ones Pa = Pa — Qi and 
Pb = Pb- q2- 

[l-Xi)b [l-X2)S 

where PaPb = S/2. In QMRK approximation one neglects the terms proportional to pb in 
p'^ and proportional to Pa in p^. Now 

Ql = XiPa + gi±, q2 = X2Pb + g2± 

and ql = qjj^ and qj = ql^ 

The explicit expression for F depends on the nature of incident particles. For example, 
if the incident particle a is a gluon the vertex has a from 

n = 2/l,we"(Pa)e"'(K), 

where /*^/ is gauge algebra structure constant and e is the gluon polarization vector. For 
quark scattering 

T\ = 2t\^,u{p^)^v{p':), 
where t\j^i is now a matrix in fundamental representation. 



In general, color structure of (|A.1| ) can be presented as j'^iy^i'z j-^^a -^j^gj^g different 
T's are color algebra generators in appropriate representation (. . . denote color indices of 
particles produced in the central rapidity region). The corresponding factor in the cross 
section reads 

tr(T^iT5)tr(r^^2T5i^2)tr(Tfe^2TfJ^) 

Using the well known property of irreducible representations , tr (T^^^Ta^) ~ the 
summation over final and averaging over initial color indices can be converted into the 
averaging over ii and 12 indices in MJ^^^^^ in (|A.1| ) with appropriate factors included into F's 
(and thus into structure function definition). As these additional factors are completely 
independent of the structure of M we are having unambiguous determination of the 
factorization of the cross section into structure functions and a generalized cross section 
for the scattering of virtual particles described by them. The correct normalization is, in 
particular, crucial for getting a correct limit of coUinear factorization in which the hard 
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cross section for the scattering of on-shell particles described by the usual ("integrated") 
structure functions should have correct color factors built in. 

When particles produced include gluons the amplitude M^^^j gets contributions not 
only from diagrams of 2 — > n + 2 type with n lines attached to the t-channel gluon, but 
also from diagrams with bremsstrahlung from pa [p'a) ^ind pb (p^) lines . These can be 
written in a form of ( |A.1| ) but, having no gluon with momentum qi (^2) in the t-chanel, 
give contribution to M proportional to qlj_ Thus the amplitude M^^^^ has the form 

Note that if even one of n particles is a gluon produced by bremsstrahlung from pa (or p'^ 
the corresponding diagram contributes to M*^^\ In the collinear factorization limit only 
M*-^^ contribution survives. 

Our next goal is to show that the amplitude M can be rewritten in such a way, 
that the nonsense polarizations dominating the fluxes coming to the hard vertex can 
effectively be traded for the transverse ones providing a basis for interpreting the hard 
block contribution as a (modified) cross section. To do this let us consider the amplitude 
742^„+2 in the axial gauge with the gauge vector lying in {pa,Pb) plane, n = apa + bpb- In 
this gauge F's do not change and the numerator of the gluon propagator is 



df,v{q) = g^w — — \-n 



V 



n ■ q [n ■ qY 

It is straightforward to check that the following important relations hold 

Padt^uiqi) = — -qi±,u, Pbdf,uiq2) = -—q2L,u- (A.2) 

Xi X2 

Let us now inspect how the structure of p^^ M^^^^p'^'^ changes in this gauge 

- P'a'd;\iqM%d-;iiq,)p^,\ 

P'a'M^^ly,- p^,^d';\iq,)Mj^lp^,\ (A.3) 



b 

-(4) „M2 Ml /Vf (4) 



where M*^*) have to be calculated in the new gauge. Now one can present M^'^\ M^^^ and 
M(^) as follows 



/vf(2) ^ p7 , .(2) 

^"MiM2 g ^"l^l/i2 

V2 

/1>(3) ^ g2±,M2 ,^(3) Pb_ 
^"AtiM2 ^"Mii'2 g 

yrv(4) ^ q2J.,^l2 p7 jyjjA) pV 

^^^^ X2 xi S "1"^ 5- 
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Using now ( [A .21 ) we obtain 

Xi X2 



where 



The amphtude A2^n+2 projected onto the physical states of incoming and outgoing 
particles is, of course, gauge invariant. While gauge transformations do not change 
(when transforming from the covariant to the axial gauge) the pJ^^M^^^^Pfe^ projected onto 
the physical polarizations of outgoing particles also remains the same (see ( |A.3|) ). This 
proves one can rewrite M^^^j in the form where the t-channel gluons having momenta 
qi± and q2± are having transverse polarizations in the original covariant gauge when the 
amplitude is projected onto the physical subspace. 

Appendix B Cross sections of pair production in high 
energy factorization 

Let us introduce the following notations 

s = 2{kik2ch.{Ay) - ki±k2±); 

t = -iqi± - ki±f - kik2e'^y, u = -{qi± - k2±f - kxk2e~^^\ 
S = X\X2S = kl + kl + 2kik2ch(Ay) , 



where ki = \/klj_, k2 = \ k^^ and ki±k2± is a dot product with 2d Euclidean metric 



1_L5 1^2 — y '^2± 

The combined contribuion from gluons and quarks (fermions) to gg scattering has the 



form (see (uA 



A. Agluons ~l~ A fermions (-^■-^) 



B.l gg gg 

•A-gluons A-\ ~l~ ^2 



A, 



1 I 1 qhl e'^y 

tu Uuklkl 4:tkik2 



-Ay 



1 



Aukik2 Akfkl 



+ 



(l + -)ch(At/)- 

s 



ki 

k2 

k2 



-#[(l + ^e-^^)U(l + ^e'^^)- 
4s k2 t ^ ki 'u 



(B.2) 
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4 ^ ih^ - + hk^e-^y s J [ (k,^ - + k.k^e^y sjj' 

E = (gi^ - g2±)(fci± - k,^) - - ql){kl - kl) + 2k,k,s\i{^y) [l - 
B.2 gg qq 

2 



A. f ermions "^1/ ~^ "^2/ 



^??2V. , , , ^A K\ fih±-qi±)\k,^-q,^f-klk'^' 



A., - {2^-^[l + k.k2sHAy)i---))-[ -- - - J + 

1 / (fc2± - gi±)^ - k,k,e-^y E\ / (fcix - gi^)^ - fcifc^e^^ E\ | 

2 V ik2± - qii_Y + k^k^e-^y s J [ (k,^ - q^^f + k^k2e^y ^ sj] ^ ' > 

and 

A2f = 



{ki± - qi±f{k2± - qi±f - k^k^V q\q\ 



tu ) tu 

where E is the same as for gluons. 

In multy-Regge kinematics only term leading in Ay oo survives and 



AmRK — AMRK,gluons — ,21,2 

Kl K2 
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Figure captions 



The cross sections depicted in figures are calculated for unintegrated structure func- 
tions CTEQ5M and AKMS |^ with constant = 0.2. Gluon and quark contribu- 
tions (with Uf = 5) are added up. 

Fig. H] Cross section Eq. (|T^) calculated for = 2 GeV, yo = and Ay e [0, 2] 
a with CTEQ5M structure functions, = 1.8 TeV 
b with CTEQ5M structure functions, y/S = 14 TeV 
c with AKMS structure functions, \/S = 1.8 TeV 
d with AKMS structure functions, = 14 TeV 

Fig. ^ Angular asymmetry p(A0) for = 2 GeV, yo = and Ay = 0.5, 1, 2. Specific 
parameters for a, b, c, d same as above. 

Fig. ^ Angular asymmetry p(A0) for 7/1^2 = ±0.5 and ko = 4, 8, 12, 16 GeV. Specific pa- 
rameters for a, b, c, d are the same as earlier. 

Fig. ^ Momentum asymmetry Eq. (p!5[) for yi 2 = ±0.5. Specific parameters for a, b, c, d 
same as above. 



16 





Figure 1 
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Figure 3 
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Figure 4 
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